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ABSTRACT 

_  The  behavior  of  coat-substrate  systems  subjected  to  repeated  impirige- 

wents  of  liquid  droplets  was  investigated.  The  systems  studied  consisted 
of  a  thick  homogeneous  substrate  covered  by  a  single  layer  of  homogeneous 
coating  of  arbitrary  thickness.  Ba«ed  on  the  uniaxial  stress  wave  model, 
the  variations  of  the  stresses  with  time  were  determined  both  in  the 
coating  and  in  the  substrate.  Employing  the  fatigue  theorems  established 
for  the  rain  erosion  of  homogeneous  materials,  algebraic  equations  were 
derived  which  describe  the  incubation  period,  and  the  mass  loss  of  the 
coating  past  the  incubation  period.,  in  terms  of  the  properties  of  £he 
droplet,  the  coating  and  the1  substrate.  The  results  were  compared  to 
available  experimental  data  and  good,  agreement  was  found  between  the 
present  analytical  results  and  the  data.. 

The  differences  between  the  yniaxial  stress  wave  and  the  uniaxial’ 
strain  wave  models  were  also  evaluated  by  calculating  according  to  both 
models  a)  the  stress  at  the  coat-liquid  interface,  b)  the  stress  that 
would  occur  in  the  substrat  .  in  the  absence  of  the  coating,  and  c)  the 

stress  in  the  coating  after  thh  first  wave  reflection  from  t’'e  substrate. 
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SECTION  I 


INTRODUCTION 

Components  of  high  speed  aircraft  and  missiles  may  experience  heavy 
damage  when  subjected  to  repeated  impingements  of  rain  droplets .  The 
damage  to  nohmetallic  components,  such  as  plastic  radomes,  may  be  parti¬ 
cularly  severe.  To  protect  such  surfaces  from  rain  erosion,  th  *se  sur¬ 
faces  are  frequently  covered  with  a  thin  layer  of  coating.  Considerable 
research  has  been  performed  in  the  past  to  select  the  most  suitable  coat¬ 
ing  material,  and  to  determine  the  behavior  of  various  coat-subatrate 
systems  undergoing  liquid  impingement. 

The  majority  of  the  previous  studies  of  rain  erosion  of  coated  mat¬ 
erials  have  been  experimental  in  nature,  with  the  bulk  of  prior  research 
concentrating  on  the  measurement  of  an  erosion  parameter  (e.g.  weight 
loss)  under  specific  conditions  (References  1-6) .  These  experimental 
studies  provide  information  on  the  behavior  of  a  given  coat-sub3trate 
combination  under'  a  given  condition,  but  fail  to  describe  material  be¬ 
havior  beyond  the  range  of  the.  experiments  in  which  they  were  obtained. 

For  the  selection  of  the  proper  materials  and  for  the  design  of  the  ap¬ 
propriate  structures  an  analytical  or  semiempirical  model  would  be  needed, 
which  would  describe  the  response  of  coat-substrate  systems  in  terms 
of  the  relevant  parameters.  These  parameters  should  Include  the  proper¬ 
ties  of  the  coating  and  the  substrate,  the  thickness  of  the  coating,  and 
the  impact  velocity  and  size  of  the  droplet.  In  recent  years,  'regress 
towards  this  goal  has  been  made  by  Morris  (Reference  7) ,  Engel  ,  \d 
Piekutowski  (Reference  8)  and  by  Conn  and  his  coworkers  (References  9-11), 
who  analysed  the  stress  history  in  various  coat-substrate  systems.  A1- 
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though  the  results  of  these  investigations  further  our  understanding  of 
the  processes  which  contribute  to  the  failure  of  the  coating  and  the  sub¬ 
strate,  as  yet  they  are  not  capable  of  correlating  fully  the  existing 
data  and  generalizing  the  results  obtained  from  a  few  experiments. 

The  objective  of  this  investigation  is  to  develop  a  node!  which  is 
consistent  with  experimental  observation  and  which  predicts  quantitatively 
"erosion”  of  coated  materials .under  previously  untested  conditions.  In 
particular,  the  model  proposed:  here  Is  aimed  at  describing  a)  the  "incu¬ 
bation  period",  i.e.  the  time  elapsed  before  the  bims  loss  of  the  coating 
becon»3  appreciable,  arid'  b)  the  degradation  of  the  coating  past  the 
incubation  period,  as  manifested  by  its  mass  loss.  The  model  is  based  oh 
fatigue  concepts  (e.g.  References  12,  13),  and  la  along  the  lines  devel¬ 
oped  previously  for  homogeneous  (uacoated)  materials  (Reference  13). 

The  success  of  thi3  model  in  describing  the  damage  of  homogeneous 
materials  warranted  its  extension  to  coated  materials. 
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SECTION  II 


THE  PROBLEM  , 

The  problem  investigated  is  the  following.  Spherical  liqi&d  drop¬ 
lets  impinge  repeatedly  upon  a  plane,  semi-infinite  material  c-iisisting 
of  a  homogeneous  substrate  covered  by  a  homogeneous  coating  (Fig.  1). 

v 

The  thickness  of  the  coating  is  h.  The  substrate  is  taken  to  le  semi¬ 
infinite  normal  to  the  plane  of  the  surface  (x  direction  in  Fig.  1).  The 
coating  and  the  substrate  are  characterized  by  the  following  properties: 
density  p,  speed  of  sound  C,  modulus  of  elasticity  E,  Poisson’s  ratio  v, 
ultimate  tensile  strength  and  endurance  limit  Oj.  Parameter  related 
to  the  coating  and  the  substrate  are  denoted  by  c  and  s;,  respectively. 
Parameters  related  to  the  droplet  are  identified  by  the  subscript  L. 

A  perfect  bond  is  assumed  between  the  coating  and  the  substrate,  i.e. 
at  the  interface  (x=h)  the  stresses  and  the  displacements  are  the  same 
in  the  coating  and  the  substrate.  Furthermore,  the  stress  wave  propa¬ 
gating  through  the  coating  and  the  substaate  are  considered  to  be  one 
dimensional,  propagating  normal  to  this  surface  (compression  waves).  -Waves 
parallel  to  the  surface  (shear  waves)  are  neglected. 

The  diameter  of  the  droplets  d,  the  angle  of  incidence  0,  and  the 
velocity  of  impact  V  are  taken  to  be  constant.  The  spatial  distribution 
of  the  droplets  is  considered  to  be  uniform.  Accordingly,  the  number  of 
droplets  impinging  on  unit  area  in  time  t  is  (Reference  13) 

n  **  (Vcos0)qt  (1) 

where  q  ia  the  number  of  droplets  per  unit  volume.  Rain,  falling  with 
constant  terminal  velocity  is  usually  characterized  by  a  parameter  I 
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called  "intensity"  (with  units  o£  length/time)  which  is  related  to  q  By 
the  expression 


6  I 

q  r  Vi3' 


(2) 


Equations  (1)  and  (2)  may  be  combined  to  yield 


h  ,  1  (V  ccs 9) I 

*  vtd3 


(3) 


The  impingement  rate  is  assumed  to  be  sufficiently  low  sc  that  all  the 
effects  produced  by  the  impact  of  one  droplet  diminish  before  t!-e  impact 
of  the  next  droplet  (References  13,  14).. 

The  pressure  within  the  droplet  varies  both  with  position  and  with 
time.  For  simplicity.,,  the  pressure  at  the  liquid-surface  interface  is 
taken  to  be  constant.,  its  value  being  given  by  the  water  hammer  pressure 
(Reference  15) 


V  cos0 


1  + 


hCL 

pcCc 


(4) 


Although  more  accurate  representation  of  the  pressure  is  possible  (Refer¬ 
ence  15)  the  accuracies  afforded  by  the  use  of  equation  (4)  will  suffice 
in  the  present  analysis*  The  duration  of  the  pressure  at  the  interface 
is  approximated  by 


h 


(5) 


The  forces,  created  by  the  repeated  droplet  impacts,  damage  the  material 
as  manifested  by  the  formation  of  pits  and  cracks  on  the  surface,  and  by 
weight  loss  of  the  coating  material.  Experimental  evidence  indicates 
that  under  a  wide  range  of  conditions  the  woight  loss  W  varies  uith  time  t 
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as  shown,  schematically ,  in  Fig.  2a.  For  some  period  of  time,  referred 
to  as  incubation  period,  the  weight  lo$s  is  insignificant.  Between  the 
end  of  the  IhciAation  period  t^  and  a  tise  denoted  by  the  weight  loss 
varies  nearly  linearly  with  time..  After  the  relationship  between  W 
and  t  becomes  more  con? lex.  Here,  we  will  be  concerned,  only  with  the  be¬ 
havior  of  the  material  up  to  time  t^.  In  most  practical  situations  the 
usefslness  of  the  material  does  not  extend  beyond  t^. 

It  is  advantageous  to  replace  the  total  weight  loss  of  the  sample 
by  the  mass  loss  per  unit  area  m,  and  the  time  by  the  number  of  droplets 
impinging  upon  unit  area  n.  In  terms  of  the  parameters  m  and  n,  sche¬ 
matic  representation  of  tl?e  data  Is  given  in  Fig.  2b.  It  is  now  assumed 
that  the  data  can  be  approximated  by  two  straight  lines  as  shown  in 
Fig.  2b,  i.e. 

a  *  0  0  <  n^  (Sa) 

m  «  a  (n-n^)  n^<  n  <  n^  (6b) 

Thus1,  the  material  loss  m  produced  by  a  certain  number  of  impacts  n,  can 
be  calculated  once  the  incubation  period  n^  and  the  rate  of  subsequent 
mass  loss  (as  characterized  by  the  slope  a)  are  known.  Therefore,  the 
problem  at  hand  Is  to  determine  the  paramett a,  and  the  latter 
being  the  upper  limit  of  validity  of  equation  (6b).  It  is  noted  here  that 
the  above  model  is  valid  only  if  there  is  an  incubation  period.  Problems 
in  which  even  one  impact  results  in  appreciable  damage  will  not  be  con¬ 
sidered. 


4.  ster  it&iS&j&stiZfiZfe. 


lit  order  to  establish  and  a,  the  stress  history  in  the  coating 
aust  be  known.  Thus,  first  expressions  are  derived  which  describe,,  in 
suitable  fora,  the  variation  of  the  stress  with  tine  in  the  coating  and 
in  the  substrate. 


SECTION  III 


STRESS  HISTORY  OP  THE  COATING  AND  THE  SUBSTRATE 

The  variation  of  the  stress  with  tine  say  be  evaluated  by  consider¬ 
ing  either  uniaxial  stress  waves  (References  10,  11)  or  uniaxial  strain 
waves  (References  6,  7)  propagating  through  the  coating.  As  will  be 
shown  in  Section  VIII  these  two  approaches  yield  similar  results.  The 
present  calculations  are  based  on  the  uniaxial  stress  wave  model. 

When  a  liquid  droplet  impinges  upon  the  surface  of  the  coating,  a 
stress  wave  propagates  through  the  coating  (see  Fig.  3).  The  magnitude 
of  this  initial  stress  wave,  denoted  by  cr^,  is  identical  to  the  hydro-' 
static  pressure  P,.  i.e. , 


°1  *  P 


(7) 


P  is  given  by  equation  (4).  At  the  coat-substrate  interface  a  portion  of 
the  stress  wave  is  transmitted  into  the  substrate  while  a  portion  of  it 
is  reflected  back  into  the  coating,  Thus,  there  is  a  "left"  traveling 
wave  in  the  coating  of  magnitude  (M-g*  3) 

o2  "  ai  +  (3) 


In  equation  (8)  a ”  represents  the  magnitude  of  the  reflected  wave  which 
may  be  expressed  as  (Reference  8) 


Z  -Z 
s  c 

*1  Z  fZ 
s  c 


(9) 
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In  the  tine  interval  t»C  /2h  the  "left"  traveling  o,  wave  reaches  the  coat- 

C  4 


liquid  interface  and  a  new  "right"  traveling  wave  of  magnitude  °3  is 
generated  at  the  x*0  surface 


°3  *°2  +0r 


(10) 


where  a  r  is  the  reflected  wave  from  the  surface  of  the  coating  (Refer¬ 


ence  8) 


Z,-Z 


Lt  C 


2  I'+Z 
•L  c 


(ID 


In  equations  (9)  and  (11)  Z  is  the  impedance  of  the  material 


Z  =  pC 


(12) 


Introducing  the  notation 


fee 


Z  -Z„ 
.,s  c 


Z  +Z 
s  c 


(13) 


ZT-Z 
L  c 


Lc  "  Z.+Z 
L  c 


(14) 


the  magnitudes  of  the  "left"  and  "right"  traveling  waves  become 


°i ' p 


a2  *  °1  +  °1  aia+*ac> 


(15) 


0,  *  <7  (1  +  *  +  tr  J 

J  1  SC  SC  Lc 


aA  "  +  *sc  +  K,  *T.„  +  ^To 


4  ”1 

etc. 


sc  TLc 


rsc  Lc  sc' 
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Equations  (15)  may  readily  be  generalized  to  the  following,  fonts 


°2k 


1  +  * 


sc 


1  i>T 
sc  Lc 


t1  -  «.Ac>  > 


(16) 


°2i-i  .  !?k  .  ^  ,  f-i 

au  d,  sc  sc  Lc 


(17) 


'1  “1 

where  lc  is  an  integer,  k  -  1,  2*  3,  ....  . 

Note  that  the  stress  history  in  the  coating  depends  oh  the  relative 
magnitudes  of  Z  ,  Z  and  Z  «  This  is  illustrated  in  Pig.  4,  where  the 

L  C  3 

variation  of  the  stress  with  tine  is  shown  for  the  four  possible  combina¬ 
tions  of  Impedances^  After  a  long  period  of  time  (i.e.  after  a  large 
number  of  reflections,  k-n»)  the  stress  at  both  on  the  surface  of  the 
coating  (x-0)  and  at  the  coat-substrate  interface  (x-h)  approaches  the 
constant  value 

.1  +  t 

O-ju  "  1 

k ■** 


a»0  lim  o  - 
1  _  2k 


sc 


1  +  Z.  U 

.  u  C 


r.  i  +  zT  n 

sc  TLc  L  s 


(18) 


ow  is  the  stress  that  would  occur  in  the  substrate  if  the  droplet  would 
impinge  upon  it  directly  in  the  absence  of  a  coating  (see  Appendix  I) . 

It  is  evident  from  Fig,  4  that  the  coating  reduces  the  stresses  in  the 
substrate  only  if  the  appropriate  coating  material  (i.e.  appropriate  cpmr 
bination  of  ZT ,  Z  and' Z  )  is  selected  (Figs.  4c  and  4d)..  For  certain 
combinations  of  coating  and  substrate  the  mean  stresses  in  the  substrate 
are  actually  higher  with  the  coating  than  without  it  (Figs.  4a  and  4b). 
Thi3  result  clearly  indicates  the  importance  of  the  proper  selection  of 
the  material  used  as  coating  for  a  particular  substrate. 
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Variation 


Equations  (16)  and  (17)  describe  the  variation  of  the  stress  with 
time  in  the, coating.  For  our  further  calculations  it  is  convenient  to 
replace  the  stepwise  variation  of  the  stress  by  a  continuous  function. 
To  accomplish  this,  equation  (16)  is  rewritten  in  the  fora 
a 


°2k  "2  k-1 

°1  al  &i  °l  Vsc'Lc 

Equation  (19)  is  now  approximated  by  the  expression 


(19) 


a 

2k  » 


CI  °1 


( 


o  o 

«o 


2,  ,  <  k-1. 

o7>  <-  — > 


(20) 


1  1  e 

By  replacing  equation  (19)  by  equation  (20)  we  replace,  in  effect,  the 
stepwise  stress  function  with  an  exponential  curve,  as  illustrated  in 
Fig.  5.  In  equation  (20)  kfi  is  the  number  of  reflections  required  for 
the  stress  to  reach  63.3  percent  of  a  ,  To  evaluate  k  we  introduce  the 
condition  that  the  area  under  the  actual  (stepwise)  and  the  exponential 
curves  are  to  be  the  same.  This  condition  requires  that  the  following 
equality  be  satisfied 

r°»  o_  a  a 

-  T~>  (tA  J  1-  • 

k-1 


00 

l 


c  o  k~l  f —  v v  v*  «  « 

[™  -  (t-- - )  (\ji  ip  )  ]«  I  [ - ( — -  —-^)exp(“-“ — )  ]  dk 

°1  °i  °i  '8C  Lc  ^  ai  ai  ai  ke 


(21) 


Evaluating  the  summation  and  the  integral  in  equation  (21)  we  obtain 


e  1-  d/  ipT 
8C~Lc 


(22) 


Substitution  of  equations  (13)  and  (iA)lnto  equation  (22)  yields 


k  « 

e 


1  +  Z. /2 

la  8 


1  +  VZs 

1  +  VZs 


(23) 


In  the  absence  of  coating  Zg*Zc  and  ke*l,  which,  as  expected,  shows  that 
there  are  no  reflections  in  a  3smi-infinite  material. 
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The  tine  required  for  k^  number  of  reflections  to  occur  is  (see 


Fig.  3) 


.  2h 

t  *  V  — 

e  e  C 

c 


(24) 


and  the  nuafcer  of  reflections  during  this  tine  is 


It  a  t  — 

e  e  2h 


(25) 


Similarly ,  the  nutber  of  reflections  which  occur  during  the  duration  of 


the  impact  t^  (given  by  equation  5)  is 


“l 


«  t. 


c  c  , 

c  c  d 


t  2h  CL  h 


(26) 


It  is  to  be  noted  that  kg  is  independent  of  the  thickness  of  the 


coating  (see  equation  23),  while  k^  depends  on  h.  For  thick  coating 
(h/d  -♦?*)  k^  •*>  0  and  for  thin  coating  (h/d  -*0)  •*•••.  Thus,  the  ratio 

_kL 


(27) 


may  vary  between  xero  and  infinity.  It  is  conveninet  to  bridge  these 
two  limits  by  the  exponential  curve 


k  -  k£  [1  -  exp(-  j— ■)] 

e 


(28) 


or 


k  =  kg  [1  -  exp  (~y) ] 


(29) 


k  represents  the  average  number  of  reflections  in  the  coating.  The  var¬ 
iation  of  k  with  Y  is  illustrated  in  Fig.  6.  For  thikk  coating  k  becomes 


\/d  ■> 


(30) 
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SECTION  IV 
INCUBATION  PERIOD 


It  has  been  recognized  in  the  past  that  fatigue  plays  an  important 
role  in  the  eroeion  process  (References  12,  14,  16-21),  particularly  in 
the  “early"  stages  of  the  process,  corresponding  to  the  Incubation,  period. 
Applying  fatigue  concepts  to  the  problem  of  rain  erosion.  Springer  and 
Baxi  (Reference  13)  recently  established  a  semi empirical  formula  which 
describes  the  incubation  period  in  a  homogeneous  material;  Here,  Springer 
and  3axi’s-  analysis  is  extended  to  homogeneous  materials  covered  by  a 
single  layer  of  coating.  The  analysis  is  based  on  the  concept  that  fatigue 
theorems  established  for  the  torsion  and  bending  of  bars  might  be  applied, 
at  least  qualitatively ;  to.  materials  subjected  to  repeated  liquid  impinge¬ 
ment.  The  failures  of  bars  undergoing  repeated  torsion  or  bending  have 
been  found  to  follow  Miner's  rule  (Reference  22) 


fi  f7  f 

_ i.  x  -  4,  j.  *1 

N,  +Nv+  N 

11  q 


(40) 


where  f^,  f^  ,...f  represent  the  number  of  cycles  the  specimen  is  sub¬ 
jected  to  specified  overs  tress  levels  °el’  ae2  *  *  *  *"**  Nl*  N2>B’'Nq 

represent  the  life  (in  cycles)  at  these  overstress  levels,  as  given  by 
the  fatigue  (oft  versus  N)  curve,  is  a  constant. 

Let  us  now  consider  a  point  B  on  the  surface  of  the  material  as 
shown  in  Fig.  7.  Each  droplet  impinging  upon  the  surface  creates  a 
stress  at  point  B,  Assuming  that  the  force  created  by  the  droplet  at  it# 
point  of  impact  is  a  "point  force",  the  stress  at  point  B  due  to  any  one 
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Fig.  7.  Force  Distribution  on  the  Surface  of  the  Coating. 


droplet  is  (Reference  23) 


a  * 


F  <l-vc) 


(41), 


2nr 


where  p  Is  given  by  equation  (36).  Due  to  the  propagation  and  reflection 
of  the  stress  waves  in  the  coating  (as  discussed  in  the  previous  section) 
the  steess  in  the  coating  does  not  remain  constant,  but  fluctuates,  as 
illustrated  in  Fig,  8.  Fatigue  iife  of  the  material  is  generally  calcu¬ 
lated;  using  an  "equivalent  dynamic  stress"  (Reference  24) 


1 


>.  j 

) 


o  a  . 
a  ns 
1.  *  — —r — 

e  o  -a 
o  m 


(42) 


where  is  the  ultimate  tensile  strength  of  the  material.  In  the  pre¬ 
sent  case  o^  may  be  separated  into  two  parts  a&m  a ^  +  o£.  The  first 
part,  a'&  is  due  to  oscillations  about  the  mean  o^*  a  with  amplitude  o^  . 

The  second  part  a"  is  due  to  "oscillation"  about  the  mean  o!Vb/2,  with  a 
•a  '  id 

constant  amplitude  o’Vo/ 2.  Thus,  o’  is  not  a  constant  but  varies  with 

fit  3 

time.  For  simplicity,  we  assume  that  or*  is  a  constant  with  a  value  equi- 

fit 

valent  to  the  maximum  amplitude,  i.e. 


"i*  l°2  "  °l 


(43) 


Equations  (36)  and  (43)  yield 


o'  »  o 
a 


*Sc| 


(44) 
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Fig,  8.  The  Variation  of  the  Stress  with  Time  at  the  liquid  Droplet-Coat.  Interface. 


The  equivalent  dynamic  stresses  corresponding  to  the  two  modes  of  stress 
oscillations  just  described  may  thus  be  written  as 


(45)* 


(a/ 2)  ou 


a  -a/ 
u 


2 


(46) 


The  number  of  cycles  for  which  the  material  at  point  B  is  subjected  to  a 
given  stress  between  ^  and  oe+4fe  *8  equal  to  the  number  of  impacts  on 
a  dr  wide  annulus  located  at  r  (Fig.  7).  During  the  incubation  period 
the  total  number  6f  impacts  on  the  annulus,  is 


fi  »  n.2*rdr 


For  each  single  impact  the  number  of  stress  oscillations  in  trie  coating 

is  k  (equation  29),  The  total  nuaber  of  stress  oscillations  during  f 

imp  act  is,  therefore,  kf^.  Accordingly,  Miner’s  rule  becomes 
f  k  f 

£  fer  +  -~rr->  -  a,  CA8) 

Ni  Nt  1 
i 

where  is  the  fatigue  life  for  overs  tress  levels  at  and  is  the 
fatigue  life  for  overstress  levels  at  cr", 

Since  r  varies  continuously  from  zero  to  infinity,  equations  (47) 


and  (48)  may  be  written  as 


to  _ 

,  n.  2*r 

/  V~dr  + 


k  n.  2rrr 


r  1 
j  N" 


The  first  term  on  the  left  hand  side  represents  the  stress  oscillation 


about  a  aa/2  and  the  second  term  the  oscillation  about  a  *>a.  From 
m  m 
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equation  (41)  rdr  -is 


rdr  *  - 


2.ir 


F(l-2vc) 

2a2 


da 


(50) 


da  is  determined  by  differentiating  equations  (45)  and  (46) 
da  *  [- 


'♦scl'V  ,-l  ,  . 
- — J  d°e 


t(%c)-o)‘ 


(51) 


2  (a  N2 

V  -1 

da  »  [- - -  3  da"  (52) 


Substitution  of  equations  (50-52)  into  equation  (49)  results  ir< 


-J 


d-2v  ) 

°I  2iro.P - V 

t  1  2r  (4  a2) 


a. 


k2irniF 


N' 


da  *  - 
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i!>  (l-2v  ) 

Ysc  c 

.  .<  2 
4ir  a 

e 

_ 


N 


da" 

e 


(53) 


tt 

The  lower  and  upper  limits  of  the  integrals  have  been  changed  to  the 
ultimate  tensile  strength  ou  and  the  endurance  limit  a.^,  respectively. 

In  order  to  perform  the  integration  the  fatigue  life  N  must  be  known  as 
a  function  of  the  stress  o„.  For  most  materials  the  fatigue  curve  between 
au  and  a^.  may  be  approximated  by  (Fig,  9  ) 


N 


1  U 

Ve 


(54) 


where  b^  and  b  are  constants.  Equation  (54)  must  satisfy  the  conditions 


•  1 

for 

*e  “  au 

(55a) 

^2 

N  »  10 

for 

a  -  tfj 

(55b) 
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In  equation  (55b):,  10  corresponds  to  the  ’!knee"  in  the  fatigue  curve 
(Fig.  9).  Equations.  (54)  aid  (55)  yield 


N  «  (a  /o  ) 
u  e 


(56) 


- 

°uC  i  ' 
lo810*cfc  J  C 
xc 


(57) 


Substituting  equations  (56)  and  (36)  into  equation  (53)  and  integrating 
we  obtain 


ird  -o  . 

~4~  V  (l-vc> 


b-1  _  b-,  l 


-l)o 


(1  +  2  |*sci  k)  ■  h 


UC 


Introducing  the  definitions 


S  » 


4a  (b-1) 
u 


a  W 


4(aUc)(b-i) 

l-2v 


(5«) 


(59) 


1  +  2E  l*,cl 


(60) 


*  ird 
”i  "  ni  T 


(61) 


equation  (58)  becomes 
S 

n*  «  a. 
i  1  50 


(62) 


the  parameter  Sg  characterizes  the  "strength"  of  the  mater *al.  Thus, 
the  number  of  impacts  needed  to  initiate  damage  is  propositional  to  the 
ratio  of  the  "strength"  of  the  material  Sg  to  the  stress  a°  prc.iucsd  by 
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the  lop in gin g  droplets.  Such  a  dependence  of  n*  on  S  .  and  o°  is  reason- 

X  c 

able,  since  the  length  of  the  incubation  period  is  expected  to  increase 
with  increasing  Sg  and  with  decreasing  o°.  However,  in  view  of  the  fact 
that  equation  (62)  is  based  on  the  fatigue  properties  of  materials  in 
pure  torsion  and  bending,  one  cannot  expect  a  linear  relationship  to 
hoM  between  n*  and  Se/o°,  In  order  to  extend  the  range'  of  applicability 
of  equation  (62),  while  retaining  its  major  feature  (Namely  the  functional 

if  o 

dependence  of  n^  on  S^/a )  we  write 

S  a2  *2 

*  ,  ex  ,S  1 

n, -ai(^)  -v^r+iRpQ- 


-} 


(63) 


a  o 

where  both  a^  and  are  as  yet  undetermined  constants. 

For  a  homogeneous  material  (in  the  absence  of  coating),  the  incubation 
period  is  (Reference  13) 


*  ,S. 

-  ■«,  (?) 


(64) 


Both  P  and  o°  denote  an  average  stress  at  the  surface.  Note,  that  n* 
and  riiy  differ  only  by  the  factor  l/(l+2kj#gc| )♦  This  factor  represents 
the  damping  effect  of  the  coating. 

A  homogeneous  material  may  be  viewed  as  either  a  material  with  very 
thick  coating  (h/d  -*-«>»  k  ■*  0,  equation  30),  or  one  in  which  the  coating 
and  the  substrate  are  made  of  the  same  material  (i>  «*0,  equation  13). 

It  i3  evident  that  for  either  one  of  these  conditions  equation  (63)  re¬ 
duces  to  equation  (64),  provided  that  the  constants  a^  and  have  the 
appropriate  values,.  To  ensure  that  in  the  limits  (fc*  0  and/or  *  ■>  0) 

equations  (63)  and  (64)  become  equal  we  adopt  here  the  same  values  for 
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and  as  were  derived  by  Springer  and  Baxi  (Reference  13)  for  homo- 
*  —6 

geheous  materials.  Using  the  values  a  7, 1x10  and  afS.l  we  obtain 


-6  ' 
J.lxlO  <=f) 


(65) 


Equation  (65)  gives  the  incubation  period  of  a  single  layer  of  coating 
of  arbitrary  thickness.  The  validity  of  the  model,  must  now  be  evaluated 
by  comparing  this  result  to  experimental  data.  The  comparison  i-<*  pre¬ 
sented  in  Fig.  10.  In  this  figure  all  the  data  are  included  for  which 
both  n^  and  the  relevant  material  properties  (ay,  cr^,  b2>  v,  E,p  for 
both  the  coating  and  the  substrate  were  available.  As  con  be  se  n, 
there  is  excellent  correlation  between  the  model  and  the  data,  lending, 
support  to  the  validity  of  the  model. 

As  was  discussed  in  Section  IX,  the  present  model  is  valid  only 
when  the  incubation  time  is  greater  than  aero.  This  condition  is  met 

•k 

when  n.  >1  or,,,  according  to  equa£ion(65) ,  when  S  /o  >  8.  Thus,,  an 
1  c 

incubation  period  exists  if 
* 

ni  ” 1 

S  /o°  >  8  (66) 

e 

When  S  /o°  is  equal  to  or  less  than  8  damage  will  occur  even  upon  one 
impact  per  site.  This  is  most  likely  to  occur  at  high  impact  velocities 
in  which  case  ,o°  is  high  (since  o°  ~  P  ~  V).. 


The  value  for  the  constant  a^  was  given  in  Reference  13  as  3.  *xl0“4. 
This  value  was  obtained  by  using  the  stress  o  instead  of  oe  in  calcu¬ 
lating  the  fatigue  life.  When  a  is  replaced  by  <re  aj_  becomes 
7.1xl0“6  (see  Appendix  II), 
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Fig,  10.  Incubation  Period  nj  versus 
(Eq.  65).  Symbols  Defined 


-Present  Model 


Data; 

*  Lapp  et  al  1955 

*  Lapp  et  al  1956 
a  Lapp  et-  al  1958 

o  Schmitt  et  al  1967 

*  Schmitt  1970 


'(Jo 

Se/o”°.  Solid  Line:  Model 
Ln  Table  1. 


SECTION  V 

SATE  OF  MASS  REMOVAL 

The  bui  removal  rate  ,pf  coat-substrate  systems  can  be  calculated 
In  a  manner  analogously  to  the  mass  removal  rate  of  homogeneo. materials. 
The  analysis  relevant  to  homogeneous  materials  is  given  in  Reference  13. 
Parts  of  this  analysis  will  be  repeated  here  for  the  sake  of  o  mplete- 
ness,  and  to  enable  the  reader  to  follow  the  discussion  withoa.  the  need 
of  constant  referral  to  the  earlier  reference. 

Beyond  the  incubation  period,  erosion  of  the  surface  of  t'f*e  material 
(as  expressed  in  terms  of  mass  loss)  proceeds  at  a  nearly  constant  rate 
as  shown  in  Fig.  2b.  In  order  to  calculate  this  erosion  rate,  an  analogy 
is  drawn  again  between  the  behavior  of  the  material  upon  which  liquid 
droplets  impinge,  and  the  behavior  of  specimens  subjected  to  torsion  or 
bending  fatigue  tests.  Experimental  observations  shot  that  in  the  latter 
case  the  specimens  do  not  all  fail  at  once  at  some  "minimum  life",  but 
their  failure  is  scattered  around  a  "characteristic  life".  Fo-  specimens 
in  torsion  and  bending  tests  the  probability  that  failure  will  occur 
between  minimum  life  and  any  arbitrary  longer  life  n  may  be  estimated 
from  the  Weib.  ill  distribution  (Reference  25) 

n— n  ^ 

p  **  1  -  exp[-  (~“)  ]  (67) 

na 

where  n„  is  the  characteristic  life  corresponding  to  the  63.2  percent 

a 

failure  point  and  g  is  a  constant  (Keibull  slope).  For  (n-n^)/  «1 

equation  (67)  may  be  approximated  by 

n  _  6 
n-n. 

P  *  (-“)  (68) 

na 
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The  probability  p  can  also  be  taken  as  the  ntpber  of  specimens  that  fail 
between  h^  and  n,  If  the  material  undergoing  erosion  due  to  liquid  im¬ 
pingements  is  considered  to  be  made  up  of  many  saall  "parts”,  then  the 


aaount  of  aaterlal  eroded  (aaas  loss)  la  proportional  to  p,  i.e. 


0 


ii  *3(-n7 


n*-nj 

«,(-* — ) 

3  ni 


0 


(69) 


p  is  the  density  of  the  aaterlal  being  eroded.  In  equation  (69)  a  was 
nondlnenslonallzed  with  respect  to  p  d  In  order  to  render  the  proportion^ 
ality  constant  diaensidnless.  Equation  (6b)  is  now  rewritten  in  dim- 
enslonless  fora 


a 

pd 


»p  d3/ A 


(n*-n*) 


(70) 


Equations  (69)  and  (StO)  give 

9-1 
a 


(ri*-n*) 


(71) 


wpd3/A  3 

According  to  equation  (71)  the  mass  loss  rate  a  depends  on  the  total  num¬ 
ber  of  inpacts  n.  However,  our  model  postulates  a  constant  mass  loss 
rate  (i.e.  a  is  independent  of  n,  see  Fig.  2b),  at  least  when  h^<n<n£. 
This  requirement  can  be  met  by  setting  0*1.  Such  a  value  for  0  is  not 
unreasonable  under  high  frequency  loading  (Reference  21).  The  character¬ 
istic  life  n  is  related  to  the  minimum  life  n. ,  This  relationship  may 
a  x 


be  eppresaed  suitably  as 
*  „  a5 


na  “  Vi 


(72) 
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where  a^  end  e^  are  constants.  Introducing  the  dimensionless  uses  loss 
sate 


wpa3/4 


(73) 


equations  (71-73),  together  with  the  assumption  0-1  yield 


“  a3  — 


(74) 


The  a  given  by  equation  (74)  applies  to  both  homogeneous  materials 
and  to  coat-substrate  systems.  For  homogeneous  materials  the  ‘'Slues  of 
and  a^  were  determined  by  Springer  and  Baxl  (Reference  13)  and  were 
found  to  ba  a^-p.023  and. a^-0.7.  Similarly  as  for  the  incubation  period, 
we  adopt  the  same  values  of  these  constants  for  the  present  problem  of 
homogeneous  substrates  covered  by  a  single  layer  of  coating,  i.e. 


o  -  0.023 


(nj) 


*»0„  7 


(75) 


«,  JL  II' 

In  the  case  of  k+0  and/or  ij>  -*0  the  incubation  period  n.  reduc' s  to  n.^ 

sc  x  h 

(see  Section  IV).  Consequently,  under  these  conditions.,  o* (given  by 
equation  75)  becomes  the  same  as  given  by  Springer  and  Baxi's  formula  for 
homogeneous  materials . 

ft 

The  validity  of  the  foregoing  model  was  assessed  by  comparing  a  , 
calculated  by  equation  (75)  to  available  experimental  data.  This  compari¬ 
son,  given  in  Fig.  11,  shows  very  good  agreement  between  the  calculated 
* 

and  measured  a  values.  This  lends  further  confidence  to  the  podel. 
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Fig.  11,  Hate  of 
Solid  L: 


Dq*q: 


•  Loppetol  1955 

*  Lapp  ttol  1956 

*  Lapp.  »tal  1958 

•  Schmitt  tfoi  1967 
o  Schmitt  1970 

J _ -1— _ , _ 1_ 

io'3  I0"Z  I  o'1 

1/  n* 


osion  Versus  the  Inverse  of  the  Incubation  Period 
:  Model  (Eq:.  75).  Symbols  Defined  in  Table  1. 
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SECTION  VI 
TOTAL  MASS  LOSS  ' 

The  total  mass  loss  was  given  by  equation  (6b)  as 
m  ■  a(n-n±)  (6b) 


Introducing  the  dimensionless  parameter 


m 

P  d 
c 


equations  (6b),  (70)  and  (73)  yield 


(76) 


a*  (r.*-n*) 


(77a) 


or 


m.  _  *  * 
"  -n 

rt  A 


('7  7b) 


According  to  equation  (77b)  it  should  be  possible  to  correlate  all  ero¬ 


sion  data  on  a  m*/a*  versus  (n*-n*)  plot.  Therefore,  we  have  included 
all  the  existing  data  on  such  a  plot  (Fig.  12).  In  this  figure  the  the¬ 
oretical  result  given  by  our  model  (equation  77a)  is  also  indicated.  The 
agreement  between  the  model  and  the  data  is  quite  good,  particularly  in 
view  of  the  large  errors  inherent  in  many  of  the  measurements. 


* 
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Fig.  12. 


Coaparison  of  Present  Model  (Solid  Line,  Eq.  77b)  with 
Experimental  Results.  Symbols  Defined  in  Table  1. 
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SECTION  VII 

LIMITS  OF  APPLICABILITY  OF  MODEL 

The  results  presented  in  Sections  II-VI  are  valid  When  (a)  there 
is  a  finite  incubation  period,  and  (b)  the  mass  less  varies  linearly 
either  with  time  t  or  with  the  number  of  Impacts  n.  The  first  of  this 
condition  is  net  when  the  following  Inequality  is  satisfied  (as?  equa¬ 
tion  66) 

nA  >  1  (66a) 

According  to  equation  (65)  this  condition  nay  also  be  expressed  as 

S  /b  >  8  (66b) 

e  o 

Equations,  (66a)  or  (66b)  provide  the  lower  limit  of  the  applicability 
Of  the  model.  The  upper  limit  beyond  which  the  present  model  cannot  be 
applied  is  determined  by  the  second  condition  given  above,  namely  that 
the  mass  loss  must  vary  linearly  with  t  or  n.  An  estimate  of  this  limit 
was  made  by  observing  that  up  to  about  n*3ni  the  data  obtained  at  various 
values  of  n  did  not  show  any  systematic  deviation  from  the  model.  Thus, 
the  results  are  valid  as  long  as  the  number  of  impacts  is  less  than  three 
times  the  incubation  period,  i.e. 

n  <  3n^  (78a) 

or  in  dimensionless  form 

n  <  3n^  (78b). 
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Using  equation  (65)  we  obtain  the  following  expression  for  the 


upper  limit 


n  <21.3  x  10”6  (— ) 

5° 


Note  that  the  two  limits  expressed  by  equations  (66)  and  (78)  do  not 
inpose  any  constraints  on  either  the  naterial  or  the  iqpact  velocity. 
Thun ,  the  results  are  valid  for  any  material  and  for  Shy  velocity,  pro¬ 
vided  that  the  experimental  conditions  fall  within  the  range  specified 


by  aquations  (66)  and  ( 78) . 
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SECTION  IX  ' 

COMPARISON  BETWEEN  THE  RESULTS  OF  THE 
UNIAXIAL  STRESS  AND  STRAIN  THEORIES 

It  vai  discussed  in  Section  III  that  the  stresses  iu  the  coating  aay 
be  evaluated  "by  assuming  either  uniaxial  (one  dimensional)  stress  waves 
or  uniaxial  (one-dimensional)  strain  waves  propagating  through  the  mater¬ 
ial.  The  uniaxial  stress  wave  model  was  applied  to  the  problem  by  Conn 
et  al  (References  10,  11)  and  by  Engel  and  Plekutcwskl  (Reference  8). 

The  mi  axial  strain  model  was  employed  by  Morris  (Reference  7).  There 
has  been  considerable  speculation'  in  the  literature  (References  16.,  26;, 
21)  as- to  which  approach  yields  more  accurate  results,  Here,  we  examine 
briefly  the  differences  in the  uniaxial  stress  and  strain  models,  These 
differences  can  best  be  illustrated  using  a,  graphical  solution  method 
(Reference  7),  First  let  us  consider  the  impact  of  a  droplet  on  a  homo¬ 
geneous  (uncoated),  naterial.  Upon  impact  one  dimensional  stress  waves 
propagate  into  the  solid  and  the  liquid  with  velocities:  vg  and  v^,  re¬ 
spectively.  The  stress  at  any  point  behind  the  wave  front  in  either  the 
solid  or  in  the  liquid  is  given  by 

O  *  pvu  (80) 

where  u  is  the  particle  velocity  ct  the  point  and  p  is  the  density  of  the 
material.  The  wave  velocity  v  is  specified  by  the  relationship 

v  ■  C  -f  Bj-u  +  BjU?  _  (81) 

G  is  the  velocity  of  the  sound  in  the  material.  B.  and  Bn  are  constants. 

X  w 

The  o  versus  u  curve,  shown  in  Fig,  13,  is  called  the -Ranklne-Rugoniot 

-40-  . 

1 1*1.  <y  tf~/nrr<#  v^f~<  v.v  -aTa ■*  -v  *-*— 


for  the  coating  and  the  liquid 

d  »  p  v  u  (curve  1)  (84) 

c  c  c  o 

crj^  •  PLyL(V-uo)  (curve  2)  (85) 

are  drama  on  a  versus  uq  plot.  The  Intercept  of  these  curves  yields 
the  stress  errand  the  particle  velocity  at  surface  of  the  coating  (*■<)). 
Equation  (84)  Is  based  on  the  properties  Of  the  undisturbed  coating. 

The  Rankine  Hugonoit  relationship  for  the  coating  behind  the  stress  wave 
Is 


(a  -a.)  *  P  v  (u,-u  )  (curve  3)  (86) 

C  1  C  C  1  o 

Finally,  for  the  substrate  we  have 


a  -  pvu  (curve  4)  (87) 

s  s  s  o 

Curves  (3)  and  (4)  are  also  drawn  on  the  a  versus  uq  plot.  The  lnter-° 
cepts  of  curves  (3)  and  (4)  and  (2)  and  (4)  give  and  qmi  respectively. 
Construction  of  a  typical  a  versus  uQ  plot  is  illustrated  in  Fig.  15. 
Figure  15a  shows  the  results  for  the  uniaxial  stress  theory  (v^“CL, 
v  **C  ,  v  »C. )  for  the  oondlv !  on 

C  C  8  L 


>  o  C  <  p  C 
yc  c  Ks  s 


(88) 


For  the  tmiaxlal  strain  aodel  the  wave  velocities  v. ,  v  and  v  are  not 

L  c  s 

constants.  However,  if  the  condition 


P,vT  >  P  v  <  P  v 
L  L  c  C  S  S 
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(89) 


for  the  costing  said  the  liquid 

d  -  p jr-u-  (curve  1)  (84) 

c  c  c  o 

PL  "  PLvL(V-uo)  (curve  2)  (85) 

ere  drawn  on  a  versus  uq  plot.  The  Intercept  of  these  curves  yields 
the  stress  o^ahd  the  particle  velocity  at  surface  of  the  coating  (x»Q) . 
Equation  (84)  Is  based  on  the  properties  of  the  undisturbed  coating. 

The  Rankine  Hugonoit  relationship  for  the  coating  behind  the  stress  wave 
Is 


<0c-O  *  p  v  (“i-u  )  (curve  3) 

Finally,  for  the  substrate  we  have 


(86) 


P 

S  8  O 


(curve  4) 


(87) 


Curves  (3)  and  (4)  are  also  drawn  on  the  a  versus  uq  plot.  The  inter¬ 
cepts  of  curves  (3)  and  (4)  and  (2)  and  (4)  give  b2  and  omt.  respectively. 
Construction  of  a  typical  d  versus  uQ  plot  is  illustrated  in  Fig.  15. 
Figure  15a  shows  the  results  for  the  uniaxial  stress  theory  (VjVC^» 


v  *C  ,  v -C, )  for  the  oondiv !  on 
c  c  s  L 


y  °ccc *  p.c. 


(88) 


For  the  uniaxial  strain  no  del  the  wave  velocities  v. ,  v  arid  v  are  not 

L  c  s 

constants.  However,  if  the  condition 


P  v.  b-p  v 
I.  I.  c  c 


P  v 

S  8 


(89) 
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Fig.  15.  Impact  of  a  Droplet  on  a  Substrate  Covered  with  a  Single  Layer  of  Coating.  Calculation 
of  the  Stress  at  the  Liquid-Coating  Interface  03.,.  the  Stress  at  the  Coating-Substrate 
Interface  c?2  *  and  the  Stress  that  Would  Occur  on  the  Surface  of  the  Substrate  in  the 
Absence  of  Coating  7*..  (a)  Uniaxial  Stress  Wave  Model;  (b)  Uniaxial  Strain  Wave  Model. 
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is  satisfied  for  each  value  of  u.  then  the  Rankine-Hugonoit  curves  are  as 
shewn  in.  Fig.  15b.  Thus,  as  long  as  the  condition  ia  equation  (89)  is 
satisfied  a2  <  o^.  This -is  in  agreement  with  the  result  of  the  uniaxial 
stress  wave  model.  If  the  condition  expressed  by  equation  (89)  is  not 
satisfied  for  all  values  of  u  then  a\  may  be  larger  than  oB,  Whether 
°2  is  larger  or  smaller  than  «r®,  depends  on  the  relative  magnitudes  of 
Bl  aid  ®2  for  *he  licIuid>  the  coating  and  the  substrate.  The  conditions 
under  which  this  might  occur  cannot  be  specified  at  present  time,  be¬ 
cause  values  for  B2  and  ,B2  are  unavailable  for  most  materials. 

Plots  similar  to  those  presented  in  Fig.  15  could  also  be  drawn  for 

materials  with  different  relative  impedances  (i.e.  p  v  <  p  v  <  o  v 

L  L  c  c  s  s  ’ 

L  pcVc  PsVs  ’  PL^L  >  pcVc  >  PsVs’See  ^8*  •  However,  the  con¬ 

clusions  presented  in  the  foregoing  would  not  be  altered. 

It  is  noted  here  that  curves  (3)  and  (1)  in  Fig.  15  are  symmetfic  with 
resnect  tf  fl  *  #r  regardless  of  the  values  of  and  This  symmetry 
was  not  satisfied  by  the  Rankine-Hugonoit  plot  presented  in  Reference  7. 
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SECTION  X 
SUMMARY 


lhe  following  formulae  may  be  used  to  estimate  the  inctibaf ion  time 
and  the  mass  loss  of  the  coat  material  of  coat-substrate  sysfre  ,<  subjected 
to  repeated  impingement  of  liquid  droplets, 
a)  ^  Incubation  Period 


or 


-6  S„  5.7 


nf  »  7.1xiO“6 


9.05x10~6  Se. 

2  L-ol 

d  o 


5.7 


/no.  of  impact, 
unit  area  ' 


(90) 


(91) 


or 


q  Vcos0 a  o° 


5.7 


where 


4ou(b-l) 


d»2vc)(l  +  2  kj*  {  ] 


-r,  P7CT  VcosO  1  +  ill 
<j°  «  -L-  L  __ _ sc 


(time) 


1  -i* 


°lCL 

°cC= 


X  -  i ji  \bT 
scrLc 


1  *f  \li 

[1-  1-exp (-y). 

r3C  1  +  )p  Y  i 


sc 


and 


PC  -P  C 

m  ~  -S  C  C 

sc  PC  +P  C  * 

3  8  C  C 
CC  d 

y  -  c7  r  U-*a  t  1 

Ol  h  sc  Lc 


PTCt-P  C 

lb  a  -LL  C  C 

Lc  P.  C.+P  C 
L  L  c  c 


^  *  3Hr"T~  {  1  "  exP  f"  r  r  (1-*  t  )]  I 
scLc  Cr  h  sc  Lc'J  j 


(92) 


(93) 


(94) 


(95) 
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b)  Rate  of  Mass  Removal 


a*  -  92  [§H 
e 


,  -o 

o  -  70.6pcdJ  [§-] 
e 


.mass  loss.. 
^  impact  * 


and  c  are  defined  in  equations  (93)  and  (94). 


c)  Total  Mass  Loss 


*  '*  *  * 
m  «  a  (n  -nj.) 


m  *  a(n-n^). 


.mass  loss, 
^imit  area; 


Equations  (91),  (97)  and  (99)  yield  the  mass  loss  peir  unit  area  in  time  t 


m  «  70.6  p  d3  [—r]  {  (q  t  Vcos8)  - 


9.05xl0-6  ,Se 


c  lS 


(~)  }  .(100) 


Sg  and  o  are  defined  as  in  equations  (93)  and  (94). 

The  foregoing  results  are  subject  only  to  the  following  two  con- 
s  traints . 

a)  Incubation  time  must  be  greater  than  zero  (ti>0) ,  a  requirement  sat¬ 
isfied  by  the  condition 


>  7,96 


(101) 


*  " 


N>|U> 


b)  Total  tine  elapsed  nust  be  less  than  three  tines  the  incubation 
period,  i.e. 


(102) 


(103) 


* 
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TABLE  1.  Description  of  Data  and  Syabols  Used  in  Figures  10,  11  and  12 
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TABLE  1  (continued) 
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Material  properties  u3ed  in  obtaining  Figs.  10-12  are  from  References  (6),  (9)-(10^ 


APPENDIX  I 


DERIVATION  OF  EQUATION  (18) 


After  k  nmber  of  wave  reflections  the  .stress  at  the  coat-substrate 
interface  is  (equation  16) 

1  +  + 


“zk '  "i  t~~  ti-»ic*tc,kl 
sc  Lc 


(A.  1.1) 

After  a  large  number  of  reflections  (k+*.)  the  stress  approaches  the  Unit 


o«o  »  11m  a 
k-Ho 

Noting  that 


2k 


Z  -Z  Z-Z 

(-?——)  (■  ..L— .£)  <  i 

W  ASt  '  '7  -i.*7  '  T 


scrLc  vZ  +Z  7  %+Z 


4  I  W  M-  «  M 

s  c  L-  c 


we  obtain 


(A. 1.2) 


(A.l.  3) 


lim  (tf.  if.  )*+  0 
YscrLc 


k-v- 


Equations  (A.l.  1),.  (A. 1.2)  and  (A.,  1.4)  give 
o  a_,  1*H< 

-4  .  Xlm  -it-— -?$- 

°l  k~  °1  l4*sAc 


(A. 1.4) 


(A. 1.5) 


Using  the  notations  (13)  and  (14)  of  Section  III,  equation  (A, 1.5)  may 
be  written  as 


1  *  VZc 

a  «  a,  - — r— r — izr-  *  0 

11+  Zr/Z- 

Is  S 


zLvcose/(i  +  zL/zg) 
i  zTvcose/(i  +  z. /z  ) 

L  It  C 


(A. 1.6) 
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^36-^  - 


av.„;fe£  a 


We  »«  obeerve  chat  the  deaerator  «*.„««*  <*.!.«,  ls  ,<ua?  to 

“  th'  *ur,*ce  of  the  «~U*  tP-e^  ,.,  equationa  (4)  m  0)1. 

Thus,  am  Is 

ZjVcosS 


'  1  +  \<\ 


(A.  1.7) 


niie  la  the  streae  that  »ould  be  produced  on  the  surface  of  thr  sub- 
atrate  If  the  droplet  would  leploge  up*  lt  directly  («e  equation  4). 


,r*  Stof  >Sf  \  - 


— w  r  jV1  -  -  ^  “  -**  - 


-' toss' ~^y  r^xfc-  i 


APPENDIX  II 


THE  VALUE  OF  THE  CONSTANT  ^  FOR  HOMOGENEOUS  MATERIALS 


Springer  end  Baxi  (Reference  13)  calculated  the  incubation  period 
fro*  Miner’s  rule 


,£l  ,h 

N1  %  ’  \  1 


basing  on  the  stress  a  (equation  10  of  Reference  13) 
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Introducing  (sfee  equation  11  of  Reference  13) 


f  »•  n . 2*rdr 

X 


and,  (see  equation  16  of  Reference  13) 
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Springer  and  Baxi  obtained 
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Equation  (A. 2. 2)  and  (A. 2.3)  yield 


1  F(1t2y) 

rdr  «  -  ^ - do 
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Substitution  of  equation  (A. 2. 6)  into  equation  (A. 2. 5)  gives 
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(A.2.1) 


(A. 2. 2) 


(A.  2. 3) 


(A, 2.4) 


(A.2.5) 


(A.  2. 6) 


(A..  2.  7) 
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Evaluating  the  Integral  Springer  and  Has!  obtained 
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A  ni  "  al  P 


where*. 


2a  (b-1) 
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(lr2V)M-(j4-)  1 


2 a  (b-1) 

,a  _ 

1-.2V 


and  a  constant  was  introduced  in  Springer  and  Baxi's  work 
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Comparing  equation  (A. 2. 10)  with  data.  Springer  and  Baxi  deducei 

—  ^ 

values  of  a^  *  3.7x10  and  “  5.7,  i.e. 
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4  ni 


3.7x10'4(|) 


5,7 


We  compute  new  the  above  results  basing  the  fatigue  stress  N  on 
equivalent  dynamic  stress 


a  a 
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1  * 

e  o-o 
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Since  a  and  a  *>  ~  ,  equation  (A. 2. 12)  yields 
a  l  m  l 


oo 


e  2o  -0 
u 

The  replacement  of  o  by  0&  in  equations  (A, 2. 2),  (A, 2. 4)  gives 
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(A. 2*8) 
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(A.  2. 10) 
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(A.  2.11) 
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(A. 2. 12) 

(A. 2. 13) 

(A.2.14) 
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